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Abstract 

The perturbative treatment of subdominant oscillation and matter effect in neutrino beams/superbeams, 
propagating over long baselines and being used to look for CP violation, is studied here for a gen- 
eral matter density function varying with distance. New lowest order analytic expressions are 
given for different flavour transition and survival probabilities in a general neutrino mixing basis 
and a variable earth matter density profile. It is demonstrated that the matter effect in the muon 
neutrino (antineutrino) flavour survival probability vanishes to this order, provided the depletion, 
observed for atmospheric muon neutrinos and antineutrinos at super-Kamiokande, is strictly max- 
imal. This result is independent of the earth density profile and the distance L between the source 
and the detector. In the general variable density case we show that one cannot separate the matter 
induced asymmetry from a genuine CP effect by keeping two detectors at distances Li and L 2 from 
the source while maintaining a fixed ratio Li/E\ = L 2 /E 2 . This needs to be done numerically and 
we estimate the asymmetry generated by the earth matter effect with particular density profiles 
and some chosen parameters for very long baseline neutrino oscillation experiments. 
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1 Introduction 



Finding CP violation in the neutrino sector [1] is a tantalizing goal waiting to be attained [2] 
in forthcoming long baseline experiments with neutrino beams/superbeams and future ones at 
neutrino factories. The discovery and quantitative measurement of such an effect will not only open 
a new window for physics beyond the standard model, but may provide an insight into leptogenesis 
[3]. From solar neutrino studies and the reactor experiment by KAMLAND, we already know [4] 
the squared mass difference and the angle of mixing between concerned mass eigenstate pair z/ 12 
to be Amjj = m 2 — m\ ~ 7 x 10~ 5 eV 2 and sin 2 9 12 ~ 0.3 respectively (see [5] and references 
therein). Similarly, atmospheric neutrino studies [6] have revealed the corresponding parameters 
for the other mass eigenstate pair 1/2,3 to be |Am| 2 | = \m 2 — ml\ ~ 3 x 1CT 3 eV 2 and sin 2 #23 ~ 0.5. 
Other reactor experiments [7] have, however, shown that the third possible mixing angle 9 13 is very 
small, consistent with zero. These mysterious numbers have increased the importance of the task 
of getting a quantitative handle on CP violation likely to be associated with the MNS neutrino 
mixing matrix. This is sought [1, 2] to be achieved by measuring the difference between neutrino 
and antineutrino beams in the muonic to electronic flavour transition probability which can be 
obtained from the general expression 

APp a (L, E) = P[u a (0) - vp{L)\ - P[v a (0) - Vp{L)\ (1) 

over a large distance L, E being the beam energy and a, (3 being the flavour indices. If \Am\ 2 \L/E 
is chosen to be 0(1), the effect will be driven dominantly by — v T oscillation. 

A major practical problem, associated with the above task, is the occurrence of the matter 
effect in neutrino flavour transitions studied in long baseline experiments. These contribute to 
APp a (L) and induce a "fake CP violation". The latter needs to be filtered out, leaving only the 
genuine CP violation part. However, the matter effect is not merely a background but deserves to 
be studied in its own right. In particular, it can yield valuable information on the sign of Am§ 2 , 
on the MSW resonant enhancement of neutrino oscillations and on the mixing between neutrinos 
of the first and the third generations. It is important therefore to be able to compute the matter 
effect relevant to any given experiment looking for CP violation. Several studies [2, 9, 10, 11] 
have been conducted to this end. Specifically, we shall use the formalism of Arafune et al. [11] 
developed to treat the matter effect as well as the subdominant oscillation driven by Am^ and 9 12 
perturbatively to the lowest order relative to the dominant oscillation driven by |Am 2 2 | and #23 • 
Thus the evolution matrix element Sp a (x) =< vp{x)\v a (0) >, a, (3 being flavour indices, could be 
calculated to first order in Am 2 x and a, where 

a(x, E) = 2V2G F N e E = 7.56 x 10~ 5 eV 2 , (2) 

GeV gms/cc 

with N e being the electron density and p(x) the mass density of the earth, expressed as a function 
of the path length x of the beam. The parameter p was assumed in Ref. [11] to be spatially 
uniform 1 . This is not always a realistic assumption [9], especially for future very long baseline 
experiments where nonuniformity in earth's density profile cannot be neglected. One possible 

1 For instance, one could assume [12] a constant density p equal to the average density of the PREM profile [13]. 
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approach [14] is to write the density function p(x), and correspondingly a(x,E), as an average 
constant plus a spatially fluctuating part and expand the latter into Fourier modes, arguing that 
only the first few modes is important for foreseable experiments. There is, however, a large 
uncertainly [15, 16, 17] in the seismological knowledge of the latter; indeed, this feeds back into 
the average density parameter. Under the circumstances, it is worthwhile making — within the 
lowest order perturbative framework — general statements for an arbitrarily varying a(x, E) that 
can be checked by experimental measurements. This is our aim here. Our approach is basically 
analytic and the price to pay is to use lowest order perturbation theory. As will be pointed out 
in our section on numerical studies, a lowest order perturbative result is often invalid in very long 
baseline neutrino experiments. There are, however, sizable domains in the (L, E) plane where it 
should be reliable. At least, the importance of the effects considered is clearly brought out by our 
analytic considerations. 

In this paper we reexamine the treatment of Ref. [11] , assuming an arbitrary spatial dependence 
in p{x) and hence a(x,E), consistent with the approximation \a\ « |Ara§ 2 |. We give extended 
versions of the formulae of Arafune et al, accommodating such a nonuniform p(x) and in a general 
neutrino mixing basis. Even with such an arbitrary earth density profile, we find that the matter 

effect in the muon neutrino (antineutrino) flavour survival probability P[^](0) — > ^(L)] vanishes 
if the flavour conversion of atmospheric and z/ M , as observed in super-K, is truly maximal. 
However, the simple methods, proposed in Ref. [11] to filter out the matter effect terms, are 
specific to the uniform earth density assumption and do not extend to the variable density case in 
the lowest order of perturbation theory. We show that, in general, it is not possible to separate the 
'fake' matter induced asymmetry from the genuine CP effect by keeping two detectors at distances 
L\ and L 2 from the source while keeping the ratio L x / Ei = L 2 /E 2 = L/E fixed and taking a linear 
combination of P[i^(0) — > vp(L)] and P[P M (0) — > vp(L)\, as proposed for a constant matter density 
profile to the lowest perturbative order in subdominant oscillation and matter effects by Arafune et 
al. This separation has to be done numerically and in the last part we compute the asymmetries 
generated by the earth's matter effect for sample earth matter density profiles keeping specific 
experimental possibilities in mind. 

This paper is organized as follows. Section 2 contains an extension of the lowest order pertur- 
bative calculation of the evolution operator to the variable density case. In Section 3 we calculate 

the muon (anti)-neutrino flavour survival probability P[ ^(O) — > ^(L)]. This calculation is 

extended to the general flavour transition probability P[ ^(0) — > ^vp(L)] in Section 4. Our 
numerical studies are presented in Section 5. Section 6 contains a summary of our results. 
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2 Evolution operator formalism 

The formulae for the amplitude and probability of the transition of u a (0) to v${V) in vacuum are 2 

%(o)-^W) = <ML)K(o)>=Y,u* m e IE UfH , 



(3a) 



%(0)-#)) = <W ReiUciUpUZjUn) sin' 



AmlL 



AE 



j>i 



IE 



(3b) 



where Am?- = mf — m 2 and the neutrino mixing matrix U is defined in vacuum by \u a >= U a i\ui > 



, a and i being flavour and mass eigenstate indices respectively. More generally when matter is 
present, one can decompose the Hamiltonian as 



H = H + H\ 



(4) 



where H is the unperturbed part containing Am 2 2 while H' is the perturbation involving Am| x 
and a(x, E). 

In the case of an x-dependent H', we cannot use the replacement procedure of Ref. [18], but 
need to use the evolution operator formalism. Thus we write 



M x ) >= Spa(x)\up(0) >, 
where the operator S obeys the evolution equation 

dS 



i— =H S{x) + H'(x)S{x), 
ax 



(5) 



(6) 



with 



Ho = 



2E 



/0 









H\x) 



1 

2E 



\0 Am 2 2 

0' 
U | Am 2 x | U ] + 




a(x,E) 0' 





(7a) 
(7b) 



To the lowest order of perturbation, (6) can be solved by [11] 

S(x) ~ S°(x) + S'{x), 



(8) 



2 We shall use L to denote the baseline with respect to a measurement, x to denote any intermediate length and 
s for a dummy variable in any integration being performed upto x. 
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with 



S°(x) = e~ %xE \ 

S'(x) = e~ ixH H-i) [ X dse isH0 H'( S )e 

Jo 



isH° 



(9a) 
(9b) 



We turn now to matrix elements in the flavour basis and work in the approximations Am| x <C 
|Am| 2 | and \a(S, E)\ <C |Am| 2 |. First, it is convenient to define 



9(v) = ^Att&u. 

Then we can explicitly write 

S°(x)pa = ^ a + U p3 U: 3 [e-^ 2x ^- 

= 5 Pa - 2tUa 3 U* a3 e-^^ sin[g(x/E)] 
= A Pa (x/E) 

for the unperturbed part. The lowest order expression for the perturbed part 



(10) 



-i(x— s)Hq 



/3 7 



[H'(s)\ 



7<5 



-sHq 



6a 



■yi 



ex 

S'(x)p a = -i ds 
Jo 

can be rewritten, on using (7a), as 

S'{x) Pa = -i f dsE^expt-idiag {0,0, 5 (2(x - s) / E))}] u U, 

J 

[H'(s)] lS U 5j exp [-idiag {0, 0, g(2s/E)}] jj [/*.. (12) 
In the RHS of (12), H'(s) has two additive parts, one constant and one depending on s: 

H'(s) = H 1 + H a (s), (13) 

with 



1 



0' 



H 1 = —U\0 Am 2 , | U\ 
o 



'a(s,E) 0' 





(14a) 
(14b) 



Correspondingly, we can take 



S'(x) Pa = S 1 (x) f3a + S a (x) Pa: 



(15) 



where H l contributes to S* 1 and H a to S a . Since H 1 does not depend on s, we may write 

s\x)p a = -iu pi u; i {H 1 ) l5 u 5j u* aj . 



ds exp [-idiag {0, 0, g(2(x - s)/E)}] u exp [-*diag {0, 0, g{2s/E)}\ 



3,1 



.Am 



—i- 



21 



2E 



U P2 U* 2 = -iB Pa (x/E)Am 2 21 , 



(16) 



where we have used the identity 



u; t {h% u S3 



Am 



7* l- v,» -uj 2E 

Turning to S a (x)p on it is convenient to use the result 



A $i2 Sj2- 



U;, [H a {s)] l5 U Sj = 



a(s,E) 
2E 



The employment of (17) enables us to write S a (x)/3 a as 



i 

2E 

i 

2E 



rx 




/ ds exp 


-ig | 


h 





(2x - 2s)5 i3 + 2s5 



(17) 



'J3 



E 



a(s,E), 



u fH uru; i u v r*( x ,E), 



(18) 



with 

r?-(x, E) = 5 i3 5 j3 e - l ^ 2x l E ) £ ds a(s, E) + (1 - 6 l3 ) (1 - 5 j3 ) £ ds a(s, E) 

+ (1 - S i3 )8 j3 £ ds a(s, E)e~^/E) + _ ^ £ dg E) -ig{2{x - s)/E) 

(19) 

On using (19) in (18), we obtain 

S% a (x,E) = -^-S pi S la £dsa(S,E) 



E 



SiaUfaU^ / ds a(s, E) e -M*im-Mm\ S in[g(x/E) - g{s/E)\ 



l^iK 3 U 13 £dsa(s,E) e-^ E hm[g(s/E)] 



E 

~KzU^\ 2 e- i9 ^ R(x,E) 



= -iG% a (x,E). 
In (20) we have introduced the real function R(x,E): 

R(x, E) = e - i9{x/E) f ds a(s, E){1 - e l9{2s/E) } + e ig(xjE) f ds a(s, E){1 - e^* 2 '^)}, 
Jo Jo 

= 2 f X ds a{s, E) [cos g(x/E) - cos {g(x/E) - g(2s/E)}\ . 
Jo 



(20) 



(21) 



Note that R(x,E) becomes a function of x/E for a constant p(x), not otherwise. 

Finally, the flavour matrix element of the evolution operator Sp a = Sp a + Sp a + Sp a can be 
expressed as 

Sfrix, E) = Ap a (x/E) - iB f3a (x/E)Am 2 21 - iG%(x, E), (22) 

where Ap a (x/E),Bp a (x/E) and G^ a (x,E) are given by (11), (16) and (20) respectively. We can 
now calculate the transition probability Pp a to the lowest order in a and Am^ in the approxima- 
tions Am^ <C jAm^l, \a(s,E)\ <C (Am^-J stated already. We obtain 

Pp a = Pp a + Pjiai ( 23 ) 



with 

P; a = 2Im(A*, a G a pa ), (25) 



PL = A* 0a A Pa + 2Im(A*Bp a ), (24) 



3 The flavour survival probability z^(0) —> v^L) 

Using (11), (16), (20) and choosing a — ft — //, we have 

A^(x/E) = l-2t\U^\ 2 e-^ x/E) sm[g(x/E)}, (26a) 



B^x/E) = \U, 2 \ 2 ^ (26b) 

G^(x,E) = lUe3 ^ 3l \ -^R(x,E), (26c) 

with R(x,E) as defined in (21). In the matter free case, a(s,E) = and replacing x by L we 
obtain 

P^{L/E) = l-A\U^\ 2 sm 2 [g{L/E)] 

+ 4\U, 3 \ 4 sm 2 [g(L/E)] + ^\U, 2 \ 2 \U, 3 \ 2 sm[2g(L/ E)]Am 2 21 . (27) 

If we consider this transition to be overwhelmingly driven by a two flavour oscillation, as done in 
the super-K analysis [6], we can ignore the third RHS term in (27) to see that the depletion 1 — 
P® (L, E), i.e. the flavour transformation and hence mixing for a muonic neutrino or antineutrino, 
would be maximal for {U^l = l/y/2. Next, we consider propagation in matter and keep a and 
hence M(L). Then we have 

P^(L,E) = P°^(L/E) + P^(L,E), (28) 

where 

P^(L,E) = m ^\ 2 - 1) sm[g(L/E)\ R(L, E). (29) 
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The asymmetry AP W (L, E), defined in (1), is thus 

AP m (L,E) = ^ e3|2 ^ 3|2 (2|^ 3 | 2 -l) sm[g(L/E)} [R+(L, E) - R~(L, E)} 

= mZ? ^ ? W^\ 2 - l)sm[9(L/E)]R+(L,E). (30) 

In (30) R ± (L, E) is obtained from (21) by using a(s, E) = ±\a(s, E)\ respectively so that R~(L, E) = 
—R + (L,E). Thus we see that the matter effect contribution to the muonic flavour survival prob- 
ability Pfinix) vanishes if \U^\ = l/\/2. Moreover, even in matter and to the lowest order of 
perturbation, any nonzero asymmetry AP^, detected from the muon (anti) neutrino flavour sur- 
vival probabilities P( ^ — > ^) will signal a deviation from the condition for the strictly maximal 
mixing of atmospheric neutrinos at super-K. This can be used in future to sensitively probe any 
deviation of \U^\ from its maximal value l/y/2. 

4 General oscillation probability in matter: v a — > vp 

From (24) and (26a,b) we have 

PjL(L/E) ~ 6p a [l-4\U a3 \ 2 sm 2 [g(L/E)}] + 4| C/^ | 2 1 ^3 1 2 sin 2 

Re (U* 3 Up 3 U a2 U; 2 )sm[2g(L/E)} 



■+ 



E 

-2 Im (U* a3 U p3 U a2 U; 2 )sm 2 [g(L/E)] 

while (25) and (26a, c) lead us to the expression 

PL(L,E) ~ 2Im\A*(L/E)G a 0a (L,E) 



(31) 











1 




~E 


C/a3 


1 

+ E 


C/a3 



(32) 



with R(L,E) substituted from (21). 
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We can now discuss two distinct cases. 
Case 1: a = 

• a = 7^ e 

In this case the matter independent and matter dependent transition probabilities, cf. 

P° a (L/E) = l-A\U a3 \ 2 (l - \U a3 \ 2 ) sm 2 [g(L/E)\ + -^^\U a2 \ 2 \U a3 \ 2 sm[2g(L/E)\, (33a) 
^JA^) = ^l^3| 2 |^3| 2 (2|f/ a3 | 2 -l)sin[( 7 (L/ J E)] J R(L,E). (33b) 

• a = = e 

The expression for Pg e (L) is the same as the RHS of (33a) with a — e, but the matter 
dependent part is 

P: e (L,E) = ^\U e3 \ 2 sm 2 [g(L/E)]J o L dsa(s,E)sm[g(L/E) 
-k(s/E)/2] - ^\U e3 \ 4 sm[g(L/E)]R(L,E) 
+||f/ e3 | 4 (|f/ e3 | 2 - l)sm[g(L/E)]R(L,E). (34) 

Case 2: a 

Now the matter independent transition probability is given from (31) by 

P°p*(L/E) = A\U, 3 \ 2 \U a3 \ 2 sin 2 [g(L/E)]-^^ {/m£ sm 2 [g{L / E)]--Re^ sin[2^(L/E)]} ,(35) 
with 

i = U; 2 Up 3 U a2 U* a3 . (36) 
For the matter dependent part, (32) leads to three cases. 

• a 7^ e, 7^ e 

P$ a (L,E) = ^\U a3 \ 2 \U p3 \ 2 \U e3 \ 2 sm[g(L/E)]R(L : E). : (37) 

• a 7^ e, = e 

P e a a (L,E) = ^\U a3 \ 2 \U e3 \ 2 (2\U e3 \ 2 - l)sin[g(L/E)]R(L,E). (38) 

• a — e, 7^ e 

P| e (L,E) = i|f/ e 3| 2 |^3r(2|^3r-l)sin[^(L/E)] J R(L,E). (39) 
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Both the CP violating part, proportional to im£ in (35), and the matter dependent part Pp a 
change sign when one goes from neutrinos and antineutrinos. However, the former is a function 
of x/E, while the latter involves both x/E and x in an unknown way for a general earth matter 
density function p(x) = (2\^2GFN e E)~ 1 a(x, E). Only for p(x) = constant, can all the transition 
probabilities become functions of L/E and the matter dependent part can be eliminated, to the 
lowest order of perturbation, by [9] taking Li(L 2 — Li)~ 1 [APj3 a (L 2 ) — APp a (Li)] L / E fi xe d, but this 
procedure is in general invalid for a spatially varying a(x, E). 



5 Numerical estimate of matter induced asymmetry 

In this section we present some numerical results for the matter contribution to the transition 
probabilities using the formalism developed in this paper. We choose realistic neutrino beam 
energies and both realistic and notional detector baselines and discuss the conditions for the 
validity of our perturbative calculation. For numerical studies we have used the Preliminary 
Reference earth Model(PREM) [13] and another earth model - akl35-F [19]. For the validity of 
perturbation theory for all intermediate values of s we must satisfy the conditions 

a(s,E)s Am 2 31 s 

AE AE 1 ' 

•t!bBL « 1; ^f£ « 1 (41) 

where < s < L. Note that conditions in eqns (41) are required for the expansion of S(x) in 
Section 2. Neutrino oscillation experiments either look for disappearance of the initial neutrino 
beam or for the appearance of a different flavour in the final neutrino beam at the detector. Both 
these effects are maximal when sin 2 (Am| x L / AE) ~ 1 corresponding to a "peak" in the transition 
probability. The Super-Kamiokande atmospheric neutrino data demands a Am^ ~ 3 x 10~ 3 eV 2 
[6] while the mass squared difference associated with the solar neutrino oscillation has now been 
confirmed to be around Am^ ~ 7 x 1CT 5 eV 2 by the KamLAND experiment [5]. Thus for all 
realistic experimental scenarios where the beam energy is tuned to an oscillation maximum for 
Am^, oscillations driven by the solar scale is always subdominant compared to those driven by 
the atmospheric scale, i.e. the condition (42) and the second of conditions (41) are satisfied. For 
the validity of our approximation (42) we will therefore always confine ourselves to values of 

To check our other approximation concerning the matter potential we present in Figure 1 the 
comparison of the strength of the matter term vis-a-vis the term involving Am^. The solid line in 
the left-hand panel shows aZL/AE, where HZ is the "average" matter potential for a given neutrino 
baseline L defined as 

1 r L 

al — — I a(s, E)ds (44) 
L Jo 



10 




L in km 



Figure 1: The comparison of the matter term against the Am^ term as a function of the neutrino 
baseline and for different neutrino energies. The matter term is independent of energy and is 
shown by the black solid line. 

For each experimental baseline L, aZ is computed using the PREM model [13]. Since a(s, E) goes 
linearly with the neutrino energy E (cf. Eq.(2)), the matter potential term ~aZL/<AE is independent 
of the energy of the neutrino beam. It is only a function of the average matter density and hence 
of the baseline L. The dotted, dashed and long-dashed lines show Am^L/AE for E = 1, 5 and 
10 GeV respectively. We note that for neutrinos with E = 3 GeV the approximation given in Eq. 
(40) works very well for all baselines shown. However for neutrino energies equal to or in excess 
of 10 GeV our perturbative expansion for the matter term works only for smaller L and breaks 
down at higher baselines. 

In Table 1 we give a summary of the magnitude of AP a for some specific cases. In Figure 2 we 
show AP a as a function of the baseline length L for the — > v e transition and with E = 3 GeV. 
We have checked that very similar behaviour is exhibited by AP a for the transitions v e — > v T and 
Up — > v T . The solid and the dashed lines show the AP a obtained using our formalism for a varying 
density matter profile for the earth. The dashed line is for the PREM model [13] while the solid 
line corresponds to the akl35-F model [19] for the earth matter profile. The dotted line gives the 
corresponding values for a constant density earth with p = 3.28 gm/cc. The values of parameters 
used for generating the figure and the table are shown in the captions. Comparison of the solid 
and/or dashed lines with the constant density dotted line shows that the variation in the density 
profile can lead to a change in the oscillation probability. Even at L = 3000 km, we note about 
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Table 1: The magnitude of AP a for some specific values of E and L and for the different oscillation 
channels. Values of L in kms are shown in parantheses after AP a . Also shown is the value of the 
function P(L, E) defined in Eq.(32). We have taken Am 2 x = 3 x 10~ 3 eV 2 , Am 2 21 = 7 x 10~ 5 eV 2 
and \U^ 3 \ 2 = 0.5. We present results for three different values of |C/ e3 |. 

10% difference in AP a between the constant density case and the case corresponding to PREM 
or akl35-F. At smaller L we note a small difference between AP a corresponding to PREM and 
akl35-F, which just reflects the fact that akl35-F has more density fluctuation at smaller L than 
PREM 3 . We see that the effect of density variation begins to be significant for L > 5000 kms. A 
caveat is that out formulae are not quantitatively reliable if L and E are such that either (40) or 
the first of conditions (41) breaks down. (The second condition of (41) is always safely obeyed). 
Such is evidently not the case in most of Figure 1. Work is in progress to improve upon the first 
order perturbation theory in a(s, E) so that the dependence on these conditions is reduced. 

6 Conclusions 

With the confirmation of neutrino flavour oscillations both in the atmospheric as well as the 
solar neutrino sectors, the focus now has shifted to the precise determination of the oscillation 
parameters involved. This will be possible in the currently planned and future long and very long 

3 It has been pointed out in [17] that even the so called realistic earth models like PREM and akl35-F neglect 
the local density fluctuations which can have an impact on the final oscillation probability and hence on the CP sen- 
sitivity of a specific experiment. However our expressions for the oscillation probabilities in varying density matter 
are completely general and can be applied to any earth matter density profile, local variations notwithstanding. 
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Figure 2: The CP asymmetry induced by matter effect (AP a ) as a function of the baseline L 
for the Pfj, — ► v e oscillation channel with E — 3 GeV. The dashed line and the solid line give the 
AP a calculated using the PREM model and the akl35-F model respectively, while the dotted 
line is for a constant density earth with p = 3.28 gm/cc. We have taken Am^ = 3 x 10~ 3 eV 2 , 
Am 2 21 = 7 x 1(T 5 eV 2 , \U^\ 2 = 0.5 and \U e3 \ 2 = 0.01. 

baseline experiments involving conventional superbeams and neutrino factories. The determination 
of CP violation in the lepton sector and measurement of the CP phase will be the most interesting 
as well as the most challenging goal of these experiments. The matter effect induces a "fake" CP 
asymmetry even if there is no intrinsic CP violation in the neutrino sector. Thus a knowledge of 
matter effect is extremely important in all neutrino CP violation studies. The matter effect also 
helps in ascertaining the sign of the atmospheric neutrino mass squared difference Am^ which 
has a bearing on the neutrino mass hierarchy with profound theoretical and phenomenological 
implications. The effect of earth matter on the survival and transition probabilities were studied 
earlier and ways to disentangle the "real" CP from the "fake" CP due to matter were discussed. 
However most of these studies assumed a constant density of the earth matter. 

In this paper we have used the evolution operator formalism to derive the most general lowest 
order perturbative expressions for neutrino flavour survival and transition probabilities in varying 
density matter. We have worked in a perturbative scheme, where the oscillation driven by Arn^j 
is assumed to be much larger than those driven by Am^ and the matter potential a(s, E). We 
have shown to the lowest nontrivial order that for a maximal mixing in the — v T sector, the 
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matter effect in the survival probability P MM vanishes identically. We have made numerical checks 
for the validity of our perturbation approximation and conclude that as long as E % 10 GeV, 
our approximation holds, at least upto a baseline of L ~ 4000 km. Finally, we have compared 
the results obtained with the PREM and akl35-F density profiles with that for a constant density 
earth matter. Work is in progress to extend the range of validity of our theory to longer baselines. 
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